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Abstract:  The  effective  mass  of  slow  electrons  in  polar  crystals  is  investi- 
gated using  the  continuum  model  of  Pekar  and  Fr8hlich,  Pelzer  and 
Zienau.  The  Hamiltonian  is  studied  in  a Fock  representation  so  that 
effects  of  several  quanta  can  be  assessed.  The  calculations  show  that 
the  effective  mass  is  larger  than  that  predicted  by  the  one  quantum 
solution,  but  that  the  effective  mass  is  still  small  for  weakly  polar 
crystals.  For  excess  electrons  in  strongly  polar  crystals  many 
quantum  contributions  must  be  included  and  the  present  method  is 
inappropriate.  It  is  uncertain  whether  the  effective  mass  is  small 
in  this  case.  Holes  in  polar  crystals  appear  to  satisfy  the  conditions 
for  the  validity  of  the  theory  of  Landau  and  Pekar  which  predicts  a 
high  effective  mass. 


1 . Introduction 


The  behavior  of  excess  slow  electrons  in  polar  crystals  has  been  a 
subject  of  controversy  and  confusion  from  the  theoretical  point  of  view  for 
some  time.  Landau^  and  independently  von  Hippel^  have  proposed  the  idea 
that  an  electron  can  be  self-trapped  because  a slow  electron  displaces  the 
ions  about  it  which,  in  turn,  create  a field  that  prevents  the  electron  from 
escaping.  This  possibility  has  been  discussed  by  Mott  and  Gurney^  and  led 


1)  L. Landau,  Phys.  Zt.  Sowjet  2*  664  (1933). 

2)  A.vonHippel,  Z.  Physik  101,  680  (1936);  J.  Chem.  Phys.  8,  605  (1940). 

3)  N.F.Mott  and  R.W. Gurney,  "Electronic  Processes  in  Ionic  Crystals," 
Oxford  Press,  1st  ed.  p.86. 
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to  the  expectation  that  a slow  electron  introduced  into  the  conduction  band 

should  be  trapped  after  some  time  and  should  occupy  a discrete  level.  One 

should  observe  optical  absorption  lines  corresponding  to  the  raising  of  such 

a trapped  electron  into  excited  discrete  levels  or  into  a continuum.  In  fact, 

the  trapped  electrons  were  identified  with  the  F centers  studied  by  Pohl  and 

collaborators  until  indirect  arguments  removed  this  possibility  (p.  113  of 

Ref.  3).  The  optical  absorption  lines  should  be  found  experimentally  for  all 

polar  crystals  provided  enough  self-trapped  electrons  can  be  obtained,  i.e., 

if  enough  excess  electrons  can  be  injected  into  the  crystal,  and  provided 

other  trapping  processes  do  not  have  overwhelmingly  greater  cross  sections. 

3 4) 

Experimentally,  these  lines  do  not  seem  to  have  been  found.  * 

In  addition  to  the  optical  behavior  it  was  expected  that  slow  electrons 
(or  slow  holes)  would  show  very  low  mobilities  at  low  temperatures.  It 
was  thought  that  it  would  be  necessary  for  the  ions  to  assume  new  configura- 
tions in  order  that  the  electron  be  able  to  move  and  that  at  low  temperatures 
this  would  be  improbable.  This  would  imply  strongly  temperature-dependent 
mobilities.  From  the  experimental  point  of  view  the  very  great  effect  of 
trapping  by  defects  makes  it  difficult  to  obtain  information  about  self-trapped 
electrons  by  mobility  measurements. 

The  preceding  two  examples  indicate  that  self -trapping  would  have  seri- 
ous consequences  for  the  behavior  of  slow  electrons  in  polar  crystals.  In  ad- 
dition to  the  cases  mentioned  we  might  expect  an  influence  on  prebreakdown 
currents,  breakdown  strengths,  cross  sections  for  capture  by  defects,  tem- 
perature shift  of  bands,  motion  of  holes,  etc. 

The  present  work  was  undertaken  because  of  a serious  challenge  to  the 
ideas  outlined  above.  Frohlich,  Pelzer,  and  Zienau^  presented  a quantum- 

4)  W.H.Duerig,  Phys.  Rev.  86,  565  (1952). 

5)  H.  Frohlich,  H.  Pelzer,  S.  Zienau,  Phil.  Mag.  41,  221  (1950).  This 

work  will  subsequently  be  referred  to  as  F.P.T. 
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mechanical  treatment  o f the  interaction  of  an  electron  and  the  polar  optical 
modes  of  vibration  using  a model  which  is  a good  representation  of  the  actual 
situation  in  polar  crystals.  They  concluded  that  self-trapping  does  not  exist, 
and  that  electrons  can  move  freely  throughout  polar  crystals  with  an  effective 
mass  close  to  that  of  a free,  excess,  conduction  electron.  The  discrete  opti- 
cal levels  do  not  exist  according  to  this  treatment  and,  since  the  lattice  po- 
larization is  only  weakly  coupled  to  the  electron,  the  continuous  absorption 
by  free  **lectrons  should  be  computable  by  perturbation-type  methods. 

Somewhat  earlier,  work  by  Landau  and  Pekar*^  and  Pekar^  also  showed 
that  electrons  can  move  through  polar  crystals.  However,  their  calculations 
led  to  high  effective  masses  which  could  be  of  the  order  of  hundreds  of  elec- 
tron masses  for  the  alkali  halides,  in  flat  contradiction  to  the  F.P.Z.  con- 
tention. This  disagreement  is  particularly  interesting  since  both  sets  of  in- 
vestigators start  from  the  same  Hamiltonian  but  use  different  methods  of  solu- 
tion. In  the  Landau  and  Pekar  work  discrete  optical  transitions  cannot  occur, 
but  the  strong  coupling  with  the  lattice  should  lead  to  important  modifications 
of  the  continuous  absorption  of  light  by  the  electron. 

In  order  to  clarify  the  reasons  for  the  discrepancy  mentioned  and  to 
find  correct  values  for  the  effective  mass,  we  have  used  a method  which  con- 
tains the  F.P.Z.  treatment  as  the  first  stage  but  remains  valid  for  stronger 

interaction.  Accounts  of  this  method,  as  well  as  preliminary  estimates,  were 

8) 

presented  at  the  Tufts  meeting  of  the  American  Physical  Society  ' and  in 

9) 

Progress  Report  XI  of  the  Laboratory  for  Insulation  Research.  ' The  method 
takes  advantage  of  the  fact  that  the  lattice  quanta  form  a system  of  Bosons 

6)  L.D.  Landau  and  S.  Pekar,  Zh.  Exp.  Theor.  Phys.  419  (1948). 

7)  S. Pekar,  Zh.  Exp.  Theor.  Phys.  19,  796  (1949). 

8)  E.P. Gross,  Bull.  Amer.  Phys.  Soc.,  Tufts  Meeting,  May  24,  1952} 

Phys.  Rev.  88,  153  (1952). 

9)  Progress  Report  XI,  Laboratory  for  Insulation  Research,  May,  1952. 
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with  symmetric  wave  functions;  one  works  in  a series  of  configuration  spaces 
in  a way  first  introduced  by  Fock*^  and  used  by  Tomonaga'  ^ in  his  intermedi- 
ate coupling  meson  theory.  The  method  of  approximation  used  has  been  de- 
veloped by  Schweber^^  in  a thesis  on  the  configuration  space  treatment  of 
relativistic  field  theories.  It  is  based  on  a systematic  development  according 
to  the  number  of  quanta  present  in  the  wave  function. 

For  further  discussion  it  is  convenient  to  introduce  the  dimensionless 
quantity  i/j 


which  measures  the  strength  of  interaction  between  the  electron  and  the  optical 
polar  modes.  Here  n is  the  optical  index  of  refraction,  € is  the  static  di- 
electric constant,  m is  the  mass  of  the  electron  in  the  conduction  zone  and  w 

c 

is  the  angular  frequency  of  long  wavelength  longitudinal  optical  modes.  Values 
of  q are  somewhat  uncertain,  chiefly  because  one  does  not  have  data  on  . 
However,  for  typical  polar  crystals  values  of  Q are;  TIBr,  2.5  to  4.5;  AgCl, 

3.0  to  5.5;  KBr,  5.0  to  8.0. 

The  first  figure  refers  to  mc^i0.7  electron  masses,  the  second  to 
mc=i  2.0.  We  will  show  in  this  report  that  the  treatment  of  F.P.Z.  gives 
correct  results  for  the  energy  and  effective  mass  when  a < 0.5.  For  values 
of  a up  to  2 it  is  adequate  to  include  two  quantum  terms  in  the  wave  function  while 
for  q >3  inclusion  of  at  least  three  quantum  terms  is  needed  and  a treatment 
according  to  the  number  of  quanta  present  may  be  inappropriate.  Most  of 
the  calculations  are  carried  for  a one -dimensional  case  which  is  closely  re- 
lated to  the  real  case.  For  a < 2 a linear  dependence  of  the  lowest  energy 
state  and  the  effective  mass  on  a appears  to  fit  the  numerical  results  approxi- 

10)  V.  Fock,  Z.Physik  75.  622(1932). 

1 1)  S.  Tomonaga , Prog.  Theor.  Phys.  6(1447). 

12)  b Schweber,  thesis.  Princeton  ( 1 452);  Phys.  Rev.  84,  1259  ( i 9 5 i ) . 
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mately,  while  for  large  a it  appears  that,  if  higher  quantum  terms  were  in- 
cluded in  the  wave  function,  both  energy  and  effective  mass  would  vary  more 
rapidly  with  a.  This  is  relevant  to  the  discussion  of  other  treatments  which 
will  now  be  undertaken. 

Several  important  new  attacks  on  the  electron-lattice  problem  have  ap- 
peared in  recent  months.  We  shall  try  to  summarize  the  situation  as  it  ap- 
pears at  the  present  time.  Lee  and  Pines*  ^ have  set  up  the  problem,  in 
Fock  space,  as  we  have  done,  but  have  made  use  of  a variational  method,  in 
the  spirit  of  the  Tomonaga  intermediate  coupling  theory.**'  *^  The  trial 
wave  functions  are  products  of  functions  for  each  wave  vector  so  that  correla- 
tions between  successive  quantum  emissions  are  neglected.  They  find  for  the 
effective  mass  m*  = mc(l  ♦ a/b)  and  for  the  energy  depression  of  the  lowest 
state  E = — ahu>.  The  mean  number  of  quanta  present  is  found  to  be  a/2.  These 
results  agree  with  our  treatment  for  a < 3,  where  we  are  able  to  carry  out  the 
calculations,  so  that  good  results  are  obtained  for  energy  and  effective  mass 
in  spite  of  the  poor  wave  functions.  In  addition,  the  mathematical  simplifica- 
tion is  considerable,  making  the  method  very  useful.  Lf  one  inserts  typical 

values  for  a one  sees  that  for  almost  all  cases  m*  < 2m  . These  effective 

c 

masses  are  larger  than  those  of  F.P.Z.,  but  are  still  small  and  would  bear 
out  their  contention  at  least  as  far  as  the  effective  mass  problem  is  concerned. 
Another  type  of  treatment,  valid  for  weak  interactions  has  been  presented  by 
Tyablikov.*^’  *^  He  finds  m*  = mc/(l  — a/6)  and  E = — ahu,  agreeing  with  the 
intermediate  coupling  result  to  the  first  power  of  a.  These  results,  however, 

13)  T.  D.  Lee  and  D.  Pines,  Private  communication,  to  be  published. 

14)  Cf.  P.  T.  Matthews,  and  A.Salam,  Phys.  Rev.  8l6,  715  (1952). 

1 5)  S.  V.  Tyablikov,  Z.  Exp.  Theor.  Physics  2_L  16(1951). 

1 6)  S.  V.  Tyablikov,  Z.  Exp.  Theor.  Physics  ZZ_,  513  (1952). 


mast  break  down  for  larger  a.  Pekar,^  using  the  variational  method,  finds 
a trial  function  yielding  E = —0.1  a^hw.  This  will  be  discussed  further  in 
Section  2 but  for  the  present  we  note  that  for  large  a this  energy  is  lower 
than  that  given  by  the  intermediate  and  weak  coupling  approaches.  The  de- 
tailed numerical  calculations  of  this  report  indicate  that  the  linear  law  begins 
to  break  down  for  u > 3. 

Bogoliubov*^  and  Tyablikov,^*  have  developed  a new,  adiabatic 
method  for  treating  the  electron-lattice  interaction  which  is  valid  for  a » 1. 
The  results  are  E = —0.1  and  m#  = 27(a/b)^m^,  in  close  agreement  with 

Landau  and  Pekar^  and  providing  a quantummechanical  justification  of  the 
semiclassical  work  of  these  authors.  If  one  compares  the  energy  depressions 
of  the  intermediate  coupling  with  those  of  the  adiabatic  method  one  sees  that 
only  for  couplings  stronger  than  a — 10  docs  the  adiabatic  method  give  lower 
energies.  At  this  value  the  two  theories  predict  rrf*  s:  2.7mc  and  nf  ££.  210mc> 
respectively.^  For  these  strong  interactions  the  assumption  of  uncorrelated 
quantum  emissions  breaks  down  and  intermediate  coupling  theory  leads  to 
wrong  values  for  the  effective  mass.  Since  the  effective  masses  do  not 
match  at  a = 10,  it  is  difficult  to  give  precise  values  for  a somewhat  less  than 
10.  The  numerical  calculations  indicate  that  both  intermediate  coupling  and 
a treatment  based  on  a small  number  of  quanta  with  correct  wave  functions 
begins  to  break  down  for  a > 3,  so  that  at  a = 6 the  result  is  intermediate  be- 
tween m#  = 27mc  predicted  by  adiabatic  theory  and  rrf  — 2m£  predicted  by 
other  treatments.  Thus  for  strongly  polar  crystals  the  situation  is  still 
somewhat  uncertain.  For  the  lower  values  of  a characteristic  of  the  silver 

17)  N.  N.  Bogoliubov,  Ukr.  Mat.  Zhurnal  2^  3 (1950). 

18)  S.  V.  Tyablikov,  Zh.  Exp.  Theor.  Phys.  2_^,  377  (1951). 

t The  figure  m*  = 432  electron  masses  for  NaCl  given  in  reference  6 ap- 
pears to  be  based  on  a use  of  inc  between  2 and  3. 
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and  thallous  halides  all  the  theories  predict  m*ct  2m£  in  spite  of  the  fact 
that  the  adiabatic  approximation  is  no  longer  valid.  For  positive  holes  in 
the  "valence"  band  the  high  masses  mc  would  lead  to  large  values  of  a and 
thus  to  extremely  large  effective  masses. 

One  is  interested  in  asking  whether  any  of  the  theoretical  results  in 
question  can  be  checked  experimentally.  For  the  alkali  halides,  where  the 
coupling  is  greatest  and  the  effective  mass  highest,  it  seems  very  difficult 
to  determine  the  effects  of  lattice  polarization  either  by  optical  measurements 
on  free  carriers  or  by  mobility  measurements  since  the  large  number  of  traps 
and  defects  prevent  the  attainment  of  an  appreciable  concentration  of  electrons 
in  the  conduction  zones.  The  experimental  situation  appears  to  be  more  hope- 
ful in  less  polar  crystals  such  as  the  silver  halides  where  mobility  measure- 
ments by  Haynes  and  Shockley*^  and  Smith  et  al.^^  have  shown  that  elec- 
trons can  be  moved  across  an  entire  single  crystal  without  being  trapped. 
Measurements  of  the  absolute  value  of  the  mobility  and  its  temperature  de- 
pendence would  be  particularly  important.^  A high  effective  mass  should  be 
measurable  at  a low  temperature.  The  optical  modes  are  frozen  out  and  the 
chief  determinants  of  the  Hall  mobility  will  be  a mean  time  of  collision  de- 
termined by  scattering  by  acoustic  modes,  impurity  and  defect  scattering, 
and  the  effective  mass  determined  primarily  by  interaction  with  the  polar  modes. 
However,  if  present  theory  is  correct,  m#  ^ 2rn.  which  is  too  small  an  effect 
to  overcome  the  crudity  of  the  present  state  of  theory  as  regards  a calcula- 
tion of  the  mean  collision  time. 

19)  J.R. Haynes  and  W. Shockley,  Phys.  Rev.  8_2 , 935  (1951). 

20)  L.  P.  Smith,  "Semi-Conducting  Materials,"  ed.  H.Henisch,  London, 

Butterworth,  Ltd.,  1951,  p.  114. 

t Hofstadter^ ^ gives  values  at  77°K  for  AgCl,  1570  cm^/^olt  sec.; 

AgBr,  296  cm^/volt  sec.;  TIBr,  880  cm^/volt  sec. 

21)  R.  Hofstadter.  Nucleonics  4,  9 (1949). 
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19) 

The  experiments  of  Haynes  and  Shockley  ' on  AgCl  and  Engelhard  on 

Cu  refer  to  higher  temperatures  where  scattering  by  the  optical  modes 

23) 

dominates.  These  experiments  have  been  analyzed  dsing  the  Fr*ohlich-Mott  ' 
theory  which  is  a perturbation  treatment.  This  means  that  the  treatment  does 
not  even  correspond  in  accuracy  to  the  F.P.  Z.  calculation  of  the  effective 
mass;  thus  a better  method  is  needed.  It  is  possible  to  use  the  procedure  of 
this  report  to  treat  this  case,  but  it  will  not  be  done  here. 

2.  Earlier  Treatments 

Landau's  note  is  concerned  with  the  behavior  of  an  electron  in  a potenti- 
al U.  If  mUa^/h?  ^ 1,  where  a is  the  approximate  extent,  there  may  be  bound 
states.  Landau  conjectures  that  there  are  two  possible  types  of  states  of  the 
electron-lattice  system.  Ln  one  the  electron  moves  in  a relatively  undis- 
turbed lattice,  in  the  second  one  has  a localized  electron  surrounded  by  a 
strong  lattice  polarization.  If  the  second  is  energetically  preferred,  the 
claim  is  that  the  electron  must  overcome  an  activation  energy  before  it  can  be 
trapped  into  the  localized  state  so  that  at  very  low  temperatures  excess  elec- 
trons introduced  in  the  crystal  cannot  be  trapped.  These  considerations  are 
not  restricted  to  polar  crystals  but  are  subject  to  the  criticism  that  the  nature 
and  functional  form  of  the  potential  U is  not  discussed,  the  lattice  polariza- 
tion is  not  treated  as  a dynamical  system,  and  the  statement  concerning  the 
activation  energy  is  not  proved. 

Von  Hippel^  focusses  attention  on  polar  crystals  and  discusses  the 
dynamics  of  the  capture  process.  With  the  Einstein  assumption  of  a con- 
stant frequency  for  all  the  modes  of  the  optical  branch,  one  can  treat  the 
separate  unit  cells  as  independent  and  use  the  Bohr  theory  of  stopping  power, 

22)  E.  Engelhard,  Ann.  Physik  J_7,  501  (1933). 

23)  H.  Frohlich  and  N.F.  Mott,  Proc.  Roy.  Soc.  A171,  496  ^ 1939) . 
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24) 

as  has  been  done  by  Seeger  and  Teller  ' in  their  theory  of  dielectric  break- 
down. A slow  electron  of  kinetic  energy  l/lO  ev  has  a velocity  2x10^  cm/sec. 
and  a Debroglie  wavelength  \ = h/mv  = 4x10  ^ cm  and  produces  a coulomb 
field  outside  this  region.  It  is  necessary  to  introduce  the  critical  distance 
v/a,  which  for  the  case  in  question  is  2x10/6x10  '^1/3x10  cm.  Ions 
more  distant  than  this  respond  adiabatically  to  the  electron  motion  and  form 
a potential  well  for  it.  Ions  lying  between  h/mv  and  v/a  are  shock  excited 
and  receive  energy  from  the  electron.  Thus  the  particle  will  slow  down, 
decreasing  the  distance  v/a  and  increasing  h/mv.  This  deepens  the  potential 
well  and  results  in  an  effective  mass  because  of  the  changing  polarization 
carried  with  the  electron.  However,  according  to  von  Hippel,  the  potential 
well  causes  the  kinetic  energy  of  the  electron  to  increase,  so  that  v4"  increases, 
shortening  the  wavelength  and  allowing  a deeper  well  to  form.  The  stationary 
state  is  pictured  as  a localized  electron  with  a shrunken  Debroglie  length. 

This  picture  contains  substantial  elements  of  truth,  but  is  not  quantitative 
and  neglects  the  zero  point  motions  of  the  ions  and  the  possibility  of  trans- 
lational motion  of  the  "shrunken"  electron.  F.P.  Z.,  in  their  critique  of 
the  self-trapping  concept,  assume  that  the  shrinking  process  ceases  when 
Ti/mv  = \/v  - p.  They  then  estimate  the  energy  of  interaction  as  —l/i 
(l/n?  — 1 /« ) e“ta/v  2*  — 1 fl  ( 1 fc\  — l/«  ) gf(2mu»/fo)  Oi  — o tkj.  In  addition,  the 
energy  of  interaction,  according  to  this  argument  decreases  only  slowly 
with  velocity,  so  that  the  effective  mass  is  small.  F.P.  Z.'s  argument  can 
be  criticized  as  being  independent  of  the  strength  of  interaction  and  in  fact 
ceases  to  be  valid  for  very  strong  interactions,  i.e.,  the  "shrinking"  effect 
must  be  taken  into  account.  The  limitations  of  this  argument  have  their 
counterpart  in  the  mathematical  technique  used  to  treat  the  Hamiltonian  for 
the  system,  i.e.,  F.P.Z.  restrict  themselves  to  a wave  function  which  con- 


24)  R.  J.  Seeger  and  E.  Teller,  Phys.  Rev.  54^,  515  (1938);  56,  352  (1939). 
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tains  only  zero  or  one  quantum  terms. 

Mott  and  Gurney^  assert  that  self  trapping  is  possible  in  all  polar  crys- 
tals and  that  a series  of  discrete  levels  exist  because  the  potential  falls  off 
as  l/r.  This  again  neglects  the  zero  point  motions  and  quantum  fluctuations 
of  the  lattice  vibrations. 

The  preceding  discussions  are  only  qualitative.  To  give  an  adequate 
discussion  of  the  problem  one  must  treat  the  electron  and  the  lattice  particles 
as  dynamical  systems  and  write  a Hamiltonian  describing  the  particle  and  lat- 
tice motions  and  their  interactions.  This  has  been  done  for  polar  crystals  by 
Frohlich^’  in  connection  with  the  theory  of  dielectric  breakdown  and  dis- 
cussed and  improved  by  Frohlich  and  Seitz^)  and  Callen^^.  The  total  Hamil- 
tonian may  be  written  as 

^ ” ^particle  + ^interaction  Afield. 

The  particle  portion  is  taken  as  p" /Zm  corresponding  to  a particle  of  ef- 
fective mass  m^.  The  field  portion  is  a sum  over  the  kinetic  energies  of  the 
ions  and  the  mutual  potential  energies  corresponding  to  a given  set  of  dis- 
placements, the  force  constants  depending  on  the  details  of  the  interactions. 

For  small  vibrations  the  system  is  redescribed  in  terms  of  normal  modes  in 
which  the  frequencies  and  relative  amplitudes  of  the  ions  depend  in  a compli- 
cated way  on  the  interactions.  The  interaction  energy  of  electron  and  lattice 
oscillations  also  depends  in  general  on  the  detailed  distortions  occurring  during 
the  ionic  motions.  It  is  possible,  however,  to  present  explicit  expressions  if 
one  restricts  oneself  to  the  modes  with  wavelength  large  compared  to  the  lat- 

25)  H. Frohlich,  Proc.  Roy.  Soc.  A 160,  230  (1937). 

26)  H.  Frohlich  and  F. Seitz,  Phys.  Rev.  _79,  526  (1950). 

27)  H.Callen,  Phys.  Rev.  7j>,  1394  (1949). 
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tice  period.  This  will  be  a fair  approximation  for  the  case  of  slow  electrons 
in  a polar  crystal. 

For  a diatomic  crystal  there  are  three  optical  branches  and  three  acous- 
tical branches.  If  we  let  h . be  the  normal  co-ordinate  of  the  j'th  branch  and 

J 

£ the  wavevector,  we  have  for  the  energy  connected  with  the  lattice  vibrations 


H = * 2 Z <■&  ♦ -w; ^ 

field  j = l k BJ  BJ  BJ 


) • 


For  long  wavelengths  two  of  the  three  optical  or  acoustical  branches  repre- 
sent transverse  modes  and  one  a longitudinal  mode.  The  b^.  are  connected 
with  the  displacements  of  the  ions  by  the  equations 


<V*)  = 2 2 u.Kjtt 

J'l  Is 


u.(?)  = 2 2 u_. /r) 

j = i is 


where 


^+kj(r)  = ( icr^w)  bkj^+kj  {cos}  (k  * r) 

U-kj(r)  = ( XTN  ) bkj  ^ — Is  j {cTs}  (k  * r> 


and  £ are  vectors  pointing  in  the  direction  of  polarization.  They  determine 
the  relative  amplitudes  of  the  ions  and  satisfy  = and  M_  are 

the  masses  of  the  positive  and  negative  ions  and  N is  the  number  of  unit  cells 
in  the  quantization  box.  The  above  decomposition  is  into  standing  waves  and 


the  notation 
with  k. 


• ;> 


means  that  the  sin  goes  with  — k when  the  cos  goes 


In  a polar  crystal  the  perturbing  forces  acting  on  an  electron  are  of  two 
types.  There  is  first  the  distortion  of  the  electron  clouds  of  the  nearby  ions. 
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rhi •»  interaction  is  not  important  for  estimating  the  effective  n.ass  of  siov/ 
electrons  since  the  response  of  the  electron  cloud  is  rapid  and  does  not  de- 
pend appreciably  on  the  electron  speed.  Second,  the  moving  ions,  together 
with  the  electron- cloud  distortion,  produce  a long  range  dipole  component. 

This  has  been  estimated2"*  ’ for  long  wavelength  as  follows:  The  polarization 
produced  by  a normal  mode  k is  the  dipole  moment  per  unit  volume 

hfr)  - i (iW;>  - 0-w‘T»  • 151 

where  e*  differs  from  the  electronic  charge  because  of  the  dipole  contributions 
of  the  electronic  clouds.  For  long  wavelengths  the  polarization  varies  slowly 
from  one  unit  cell  to  the  next.  The  potential  produced  is  determined  from 
Poisson's  equation 

V2»kj  = 4"V  • pkj  • w 


From  this  equation  one  concludes  that  the  transverse  modes  produce  no 
perturbing  potential.  In  addition,  one  can  show  that  the  amplitudes  £ 
and  £ ^ for  acoustic  waves  of  long  length  are  such  that  is  very  small. 

For  short  acoustic  waves  the  dipole  component  is  appreciable.  One  is  thus 
left  with  the  longitudinal  optical  modes  of  long  wavelength  as  the  principal  source 
of  the  perturbing  potential  for  slow  electrons.  For  longitudinal  waves  of  in- 
finite length  we  have 


M M 

+■ 


1 v‘A  1 


~T7? 


where  p is  the  reduced  mass  of  the  ions.  The  total  perturbing  potential  is 
then  taken  to  be 


i A 


- ? V ?7?  U)  • 


(7) 
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and  the  total  Hamiltonian  (or  the  particle  and  lattice  vibrations  is 

2 , t • 2 


*7  2 {vn; 


(8) 


The  preceding  derivation  depends  very  essentially  on  the  assumption  that  the 
wavelengths  of  importance  in  the  problem  are  of  long  wavelength  so  that  sums 
over  the  shortwave  portions  give  small  contributions.  Additional  assumptions 
are  made  in  connecting  the  constants  e*  and  = with  experiment.  It  has 
been  pointed  out  by  Frohlich  and  Mott^^^  and  Lyddane,  Sachs  and  Teller^®) 
that  u>  for  the  longitudinal  optical  modes  of  infinite  wavelength  is  higher  than 
that  for  the  corresponding  transverse  modes  (reststrahlen)  because  of  the  ad- 
ditional restoring  force  due  to  the  macroscopic  field  created  by  the  motion  of 
the  ions  in  a longitudinal  mode.  Under  the  assumption  that  the  field  acting  on 
the  ions  is  E,  rather  than  an  internal  field  more  appropriate  to  the  face- 
centered  cubic  lattice,  it  is  possible  to  derive  the  equation  ^r,  where 

is  the  reststrahl  frequency.  The  effective  charge  e#  can  also  be  connected^  ^ 
with  experimentally  observable  constants  by 


, „ ,2  4 

2 (1  ire  \/  n \ 

wt  '-{ 3“ X 0 

pa  ,xt  — n ' 


The  above  Hamiltonian  is  written  in  terms  of  standing  wave  co-ordinates. 
The  b,  and  b,  can  be  shown  to  be  canonical  co-ordinates  and  momenta  and  it 

!s  is 

is  possible  to  perform  a canonical  transformation  from  these  to  a set  of 


running  wave  co-ordinates.  With 


X k = ( + 

- Mm 


^2M u)2' 


-!s; 


_ i 


)j 


> (9) 


28)  R.H.  Lyddane,  R.G. Sachs  and  E.  Teller,  Phys.Kev.  59,  673  ^1941) 


where 


-14- 


The  Hamiltonian  becomes 


(10) 


^ = zEr 


4ire 

fv 


£ ^ (XR  sink  • q + cos  k • q) 


1 (m u>2x2  + 


Y 

TT 


k) 


(11) 


This  is  the  starting  point  of  most  discussions  of  slow  electron  behavior  in 
polar  crystals.  The  basic  restriction  to  long  wavelengths  becomes  in- 
creasingly poor  for  increasing  interaction  strength  as  the  polarization  sur- 
rounding the  electron  shrinks  to  atomic  dimensions.  However,  the  Hamil- 
tonian appears  adequate  to  settle  some  of  the  basic  conceptual  conflicts  in 
the  self-trapping  problem. 

The  expression  Eq.(ll)  can  be  treated  by  classical  mechanics.  One 
writes  Hamilton's  equations  of  motion 


X £ (Xk  co#*  * *“]jfc8in*  * 

k - 


d§  = ay  = j_ 

Ht  3 £ - mc 


d.Lu- 

dt 


ay 

~S~TC. 


-4ne 


fc  17 


sin  k • q + 


Mu2X, 


(12) 


ay 


4ire 

✓v 


J * 


It  is  difficult  to  solve  these  nonlinear  equations  for  a general  set  of  initial 


. .*V 


. ,v 


conditions.  The  equations  have  2(2N  + 3)  independent  integrals  of  motion  and 
a particular  solution  is  identified  by  specifying  the  values  of  the  integrals. 

It  is  possible,  however,  to  discuss  a particular  type  of  solution  in  which 
there  is  a definite  relation  between  the  energy  integral  and  the  momentum 
integrals.  This  solution  represents  an  electron  moving  with  constant  ve- 
locity and  the  oscillators  undergoing  purely  forced  motion.  By  studying  the 
relation  between  the  energy  and  total  momentum,  one  can  assign  an  "effective 
mass"  to  the  electron  arising  from  its  interaction  with  the  vibrations.  This 
type  of  solution  should  be  the  classical  analogue  of  the  quantum  solution  for 
the  lowest  states  where  no  free  quanta  are  present.  We  take 


p = Pq  = constant 


q = t + q 

m 

c 


X,  = m 
k c 


4nc  1 


ffk  k-  p 


sin  k •(—  t + q ) 
\ m o/ 


k _ 4ire  1 

Mu  /Vk  ♦ ♦ 

k*  P , 


fq  ) 

' m no  / 


The  expression  for  X.  indicates  that  the  field  amplitudes  must  be  assigned 

!S 

large  values  for  those  waves  for  which  the  particle  speed  is  close  to  the  wave 
speed.  If  we  restrict  ourselves  to  very  slow  particle  speeds,  k • pQ/  < u> 
even  for  the  shortest  waves  (kmax)  corresponding  to  lattice  dimensions,  i.e., 
there  are  no  trapped  electrons.  For  these  very  slow  particles  we  should  ex- 
pect the  effective  mass  to  be  unity  since  all  of  the  ions  are  excited  adiabatically 
and  there  is  no  change  in  polarization  with  particle  speed.  This  may  be  veri- 
fied by  writing  the  inCiVientum 


V i 
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„ ♦(mw2x*  y2| 


I -£z-\ 

VvkMw  / 


V k Mu 

° i =“J~ 


(^-J 

' m ' 


(14) 


and  the  energy 


E = 7^  + 


(—)  — 

•'l/vk'  Mu  . 1 


V'f  V K ' 

k k • p 


f—  1 1 

e*  w 2 Vvk  RG ' 
i £ f-; — 2 

(^-) 


m 


(15) 


By  eliminating  |5  we  can  find  a relation  between  E and  \ which  for  slow 
particles  is  characteristic  of  a free  electron. 

Let  as  consider  the  case  of  particles  at  rest.  Then  pQ  = k = 0 and 

= — "2  m ) — — y • The  energy  thus  depends  on  the  details  of  the  short- 

i^rvk7  Mw^ 

wave  cutoff.  Following  F.  P.  Z.  we  can  compute  the  potential  produced  at  a 
point  r by  an  electron  at  qQ.  Inserting  the  values  for  and  from  Eqs. 

( 1 3)  we  find 

^(xb  sink.?  cosk-T) 


. , 2 cos (k  • r - q ) 

-loir  e v ° 

2 2 

VMu>  k k^ 


(16) 


1 1 


This  is  the  potential  of  a point  charge  in  a dielectric  and  verifies  the  choice 


MJ 


1.1/1  1 \ 

T " Tn\  ~2~1 

n ' 


28) 

Fr’ohlich  and  Pelzer  ' have  presented  another  treatment  of  the  behavior 


of  a classical  point  charge  interacting  with  the  polar  modes.  The  electron  is 
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assumed  to  move  with  uniform  velocity  and  the  response  of  the  oscillators 

is  calculated.  The  electron  velocity  is  high  so  that  energy  loss  associated 

29) 

with  trapping  must  occur.  Fr'dhlich  and  Pelzer  ' treat  this  phenomenon 
schematically,  neglecting  electron  recoil.  In  addition,  they  compute  the 
energy  carried  along  with  the  electron  and  find  an  effective  mass  increasing 
rapidly  with  decreasing  electron  speed. 

Deviations  from  classical  behavior  are  expected  in  two  different  ways. 
First,  one  must  consider  the  wave  character  of  the  electron.  This  leads  one 
to  expect  that  at  least  for  weak  interactions  the  cutoff  will  be  given  by  the 
Debroglie  wavelength  of  the  electron  rather  than  the  lattice  spacing.  Second, 
there  are  quantum  and  zero  point  fluctuations  in  the  field  of  the  polarization 
quanta.  Classical  mechanics  is  expected  to  become  valid  in  the  limit  of  in- 
finitely great  interactions,  but  it  is  possible  that  the  two  types  of  quantum  ef- 
fects become  unimportant  at  different  values  of  the  coupling  constant.  Thus, 
one  can  imagine  a semiclassical  treatment  in  which  one  system  is  treated 
classically,  the  other  quantummechanically . The  treatments  of  Pekar^)  antj 
of  Landau  and  Pekar*^  are  of  this  type. 

Pekar  describes  the  electron  by  a wave  function  \j/  (r)  which  produces 
a dielectric  displacement 

B(0  - J\  ^ (*)|2  if  (17) 

The  displacement  polarizes  the  lattice,  and,  in  a continuum  theory, 

^>=(4-7)  IP  • <l8> 

29)  H.Frohlich  and  H.  Pelzer,  British  Elect,  and  Allied  Ind.  Res.  Assoc., 
Tech.  Rep.  LIT,  184  (1948). 

30)  S.  Pekar,  J.  Physic-  USSR  LD.  341  (1946). 
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Thc  interaction  energy  of  electron  and  polarization  is 

V = - J P • D dT  (19) 

and  the  electron  wave  function  satisfies  the  self-consistent  Schrodinger 
equation 

V vi/z  = . (20) 

c 

Pekar  finds 

(r)  =*  0.123  a3/2  ( 1 + ar  + 0.45  s2r2)  e'ar 


and  shows  that  this  represents  tlu  optical  activation  energy.  The  thermal 
"dissociation  energy"  is  equal  to  'he  total  energy  depression  and  is  shown 
to  be 

Etho;  —0.1 1 a2fiw  . (22) 

Actually,  there  are  two  corrections  which  should  be  made.  First  the 
continuum  states  are  also  depressed.  Second,  even  in  absence  of  inter- 
action, the  continuum  lies  liu>  above  the  ground  level  so  that  one  should  add 
hw  to  the  energy  depression  to  obtain  E^.  Pekar's  expression  becomes 
valid  when  the  energy  depression  is  large  compared  to  a quantum,  i.e.,  when 
0.1  a2  » 1 or  q»3.  For  weaker  coupling,  zero  point  motions  and  quantum 
fluctuations  are  important. 

In  a later  paper‘d  Pekar  derived  this  expression  on  the  basis  of  a 
completely  quantummechanical  argument.  Starting  from  the  Hamiltonian 


, . r * • 


. • - * . 


Eq.(8)  in  the  standing  wave  representation,  he  assumed  a trial  function  of  the 


product  form 


where 


<t> (* — qk)  = n 1 4>h(qk) 


“k  ■ pk 


The  forms  of  the  functions  l^(r)  and  are  chosen  so  as  to  minimize 

the  energy.  (J)^  satisfies  a Euler-Lagrange  equation  and  is  found  to  be 

*k  = Ak'xp[-7<<,k-V)i]  ' (24) 

where  is  a normalization  constant  and  q^Q  = = 

J'D(r)  | jj*  . dr  • This  is  the  equation  of  a harmonic  oscillator 

with  its  equilibrium  position  displaced  by  an  amount  depending  on  the  electronic 
wave  function  \|/  (r).  Pekar  showed  that  \|/(r)  is  the  same  function  used  in  the 
semiclassical  treatment  and  that  the  total  energy  (neglecting  zero  point  energy) 
is  again  —0.11  a^hu.  The  importance  of  this  calculation  lies  in  the  demonstra- 
tion that  for  large  a the  energy  depression  varies  much  more  rapidly  than 
linearly  with  a.  The  treatment  is  based  on  the  variational  method  so  that  the 
true  energy  must  be  even  lower. 

To  find  the  effective  mass  of  a moving  polaron  in  a semiclassical  ap- 
proximation, Landau  and  Pekar^  studied  3 special,  purely  'forced'  motion. 
They  assumed  = vj/  (r  — Vt),  D = D(r  — Vt),  p - p(r  — ^t),  and  developed 
the  expectation  value  of  the  energy  in  powers  of  the  velocity.  One  obtains 

4 y\ 

m*  = 27(a/b)  m . An  attempt  by  Pekar  ' to  put  this  treatment  on  a quantum- 
mechanical  footing  is  not  convincing  but  the  more  recent  work  of  Bogoliubov*^ 


* 
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1 8) 

and  Tyablikov  ' gives  the  same  value  as  the  first  step  in  a systematic  adi- 
abatic approximation  procedure.  The  validity  of  the  serr.iclassical  procedure 
might  be  checked  by  examining  the  higher  terms  in  the  adiabatic  procedure, 
but  this  has  not  been  done.  However,  there  does  not  appear  to  be  reason  for 
doubting  the  validity  at  very  large  a where  quantum  fluctuations  are  un- 
important . 

Let  us  now  briefly  look  at  the  problem  from  the  purely  quantum- 
mechanical  point  of  view  and  in  the  approximation  of  weak  interaction.  The 
energy  levels  in  absence  of  interaction  are  given  by 


(V-  £ n^hk)2 

E(u  - 0)  = ± + 2 <nk  + 

2m  k 


(25) 


and  is  a function  of  the  quantum  numbers  X.  and  N^.  For  a one-dimensional 
case  we  may  plot  E vs.  and  find  schematically  Fig.  1.  The  point  of  in- 

tersection of  the  discrete  state  with  the  continuum  is 


4 


2m  w 
c 


or  the  critical  distance  of  F.P.  Z. 

Upon  introducing  the  interaction  the  ground  curve  is  depressed.  The 
portion  of  the  curve  belonging  to  larger  values  of  X is  closer  to  the  higher 
state  curves  and  is  therefore  depressed  more  than  the  portion  near  X = 0. 
One  thus  obtains  a changed  curvature  and  effective  mass.  For  very  weak 
interactions  it  is  possible  to  apply  second  order  perturbation  theory  and 
one  obtains  for  the  energy  levels  (dropping  the  zero-point  oscillator  en- 


ergies) 

.2  , 

E — (ahu>) sin 


■m. 


(2mchw) 


W * 


3 1)  Cf.  E.  P.  Gross,  Phys.  Rev.  84,  818  (1951). 
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These  formulas  break  down  (or  strong  interactions  and  it  is  necessary  to  in- 
clude the  effects  of  higher  quanta.  This  is  done  in  the  next  section. 


3.  Schrodinger  Equation  in  Fock  Representation 


In  accordance  with  the  discussion  a Section  2 we  shall  take  as  the 
Hamiltonian  for  an  electron  interacting  with  the  polar  optical  modes 

..2 


* = Sg^lX^sink.;.  ^.cosI.q).|  X (mo2X^  ♦ ) • (11) 


ki\  rv 


In  the  quantum  theory  it  is  possible  to  use  the  approach  of  F.P.  Z.  which 
regards  "kf  as  operating  on  a wave  function  of  the  form  \Jf(q,  X^)  and  to  re- 
place p and  1/k  by  the  operators  — lh  d/3q  and  — lh  B/BX^.  In  this  report  we 
shall  adopt  an  alternative  approach  due  to  Fock^^  which  takes  advantage  of 
the  fact  that  the  polarization  quanta  form  a Bose-Einstein  system  with  sym- 
metric wave  functions.  Fock  has  shown  that  a Hamiltonian  of  the  above  type 
can  be  regarded  as  operating  on  a wave  functional  which  has  as  components 
functions  giving  probabilities  of  various  states  in  which  a given  number  of 
quanta  are  present.  Thus  the  wave  functional 

(J)  1 Xj(k,q),  X^^i'^*4!)****  ) (28) 

has  the  following  interpretation:  I XQ(q) Jq  is  the  probability  that  an 

electron  is  found  at  q in  a region  dq  and  no  quanta  are  present:  | Xj(kj, 

q)|^  dkjdq  is  the  probability  that  there  is  one  quantum  of  wave  vector  k^ 
in  a range  dkj  and  an  electron  at  q in  a range  dq,  etc.  The  sum  ol  the  ab- 
solute squares  of  the  components  when  integrated  over  the  configuration 
ranges  is  normalized  to  unity 

With  this  representation  one  writes  the  Schrodinger  equation 
{ff—  E)  = u as  an  infinite  chain  of  integral  equations.  It  i»  first  tun.Cii- 


We  note  that  and  c£  are  not  canonical  conjugates  but,  instead,  satisfy 
the  commutation  rules  [p^.C*J=  6^^ t . In  addition,  C*  is  the  complex  con- 
jugate of  since  and  1/k  are  real.  In  terms  of  these  co-ordinates  we 
may  rewrite  the  Hamiltonian  as 


y ■£=♦  s 


Zm  + r*  7EC,  £ 


Iiw, 

- 2 • 

Is 

The  zero  of  energy  may  be  chosen  so  that  the  contribution  of  the  zero  point 
oscillations  drops  out  and  one  obtains,  following  Fock  and  Tomonaga, 

{E-p2/2m}  *(5)  = 2 k‘  ’ 

I.  * 


{E_^r_hw}  X (^q)  = 


y*  iv; 

I c 1 X(q) 


+ zl/2(m^)  £ * (kik25)  «(k2) « 


-ik2  - q 


-2  2 


‘ E - “ 2hs-)  xf\  k2  = {g(k2)  e 2 *(kl 


* 8^)  e 


i kj  • q 


X(k2q)l  . 


t • 
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This  is  the  coupled  set  of  equations  cut  off  so  as  to  include  terms  with 
two  quanta  or  less.  In  this  set  p is  understood  to  be  the  differential  oper- 
ator — i h . Let  us  introduce 


♦ ♦ 

■i  k,  • q 


f »♦  1 h * * 

^ (kt  ?)  = e X (k1  q) 


C 


2<klk2<l>  = e 


-» (kj  + k£  • q) 


X (k,  k q) 
2 1 c 


If  one  then  assumes 

X(^X^  . 


(32) 


(kj  q)  = £ j (kj)  e1  ^ ’ q 


(33) 


£2(ki  k2q)  = ^(kj  k2)c‘  K 'q 
Eq.  (3  1)  becomes 


(e"^tht)x=£  s^iKiibif  ^i(ki} 


*i 


{■-*"-£(*♦* l)2}  ^l)=^l)(7^f 


X g(kz)  £2(kj  k2) 

k i 


(34) 


{E  - 2hw-^  (\  + kj  + k2)2|  C2(kj  k2)  = |g(iT2)  £j(kj) 

+ g(kj)  £j(k2)  J * 


Let  us  first  examine  these  equations  when  the  amplitudes  of  the  two 
quantum  terms  are  set  equal  to  zero.  One  obtains  from  Zq.(34b) 


: 


• • r 
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5i‘£i>- 


, v 


e - hw  - JL(t  + k)2 

2m 


Inserting  this  into  Eq.(34a)  we  find  the  expression  for  the  determination  of 


the  energy  of  the  lowest  states 


hV 


E - hu  — (X  + k)2 


This  is  exactly  the  expression  obtained  by  F.P.Z. 

To  find  the  wave  functions  and  energy  levels  in  the  two  quantum  ap- 

proximation,  one  inserts  (,  (k.  k ) from  Eq.(34c)  into  Eq.(34b).  We  thus 

2 

find  the  pair  of  integral  equations 


Y ) C (k) 

L_,  E - 2hu> --£-(*•  ♦ k + x)Z  I 1 

„ cm  " J 


1/2  r-n  g(*)  £ (*) 

= g(k)(iL)  x + g(k)G*b)  V : i ^ ^ 

o / , E — 2hw  — — (X  + k + x) 

1 — ’ 2m  ' 


The  method  which  has  been  adopted  to  solve  Eqs.(37)  is  to  treat  Eq.(37b) 

as  an  inhomogeneous  integral  equation  and  to  solve  for  £ (k)  by  interation, 

1 

i.e.,  by  finding  the  Liouville-Neumann  series.  The  expression  for  C (I)  i. 

1 

then  inserted  in  the  right-hand  side  of  Eq.(37a)  giving  an  implicit  expression 
for  E(X).  For  an  integral  equation  of  type 


-Jbi 


-26- 

♦(*)  = h(x)  + J K(x,  y)  $(y)  dy  (38) 

with  sufficiently  regular  kernel  K(x,  y),  the  condition"*2)  for  convergence  of 
the  Liouville -Neumann  series  is 

JJ |k(J.  ;)J2  di  dy  $ 1 . (39) 


ir(k) 


For  Eq.(37b)  we  may  find  the  condition  of  convergence  by  introducing 


E — hw  — 


z. 



hz  . T . ?.\Z 


/ > E — 2hu  _ ^L-  (X.  + k + x) 
x 


t 00 

*1 


(40) 


Then 


K(krx) 

«(^)  g(*)  h/2Mcj 


(41) 


|e  - 2hu  _ (\  + k + x)2J  [e  - hw  - ^2m^'  “ 


_h 
2 Mu 


V *£*> 

/E_  ..  2,r  2 


2hu>  - ti  (£  + k + x)‘ 
2m 


4.  Effective  Mass  for  One-Dimensional  Problem 

In  order  to  discuss  the  effects  of  two  and  three  quantum  contributions 
and  to  estimate  the  effective  masses  for  the  strong  interactions  occurring 
in  ionic  crystals,  it  is  convenient  first  to  discuss  a certain  one -dimensional 
case  where  the  integrals  occurring  are  easier  to  evaluate  than  for  three-di- 
mensional. One  can  show  that  the  one-dimensional  integrals  are  certain 
averages  over  the  integrands  of  three-dimensional  integrals  and  differ  from 
the  latter  by  factors  of  the  order  of  unity.  We  shall  compare  the  F.P.Z. 
solution  with  the  more  complete  solution  for  different  strengths  of  interaction. 
At  weak  inte  actions  ((J2ir  < 1)  we  find  that  the  two  quantum  terms  give  a small 

32)  W . Schmeidler , "Integralglei  chungen."  Akademische  Verlag,  Leipzig, 

1 Q SO  n 770 

- ■ • i 
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correction  in  energy  to  the  F.P.Z.  solution  and  that  inclusion  of  three  quan- 
tum terms  gives  an  additional  negligible  correction.  For  stronger  inter- 
actions (p^ir— 2)  the  two-quantum  solution  yields  appreciable  corrections  to 
the  energy  and  effective  mass,  and  the  three-quantum  terms  give  small  cor- 
rections. For  > 3 it  would  be  necessary  to  include  more  quanta  in  order 
to  obtain  adequate  results.  It  appears,  though,  that  the  expression 
Eq  = — (J^whu  suggested  by  intermediate  coupling  theory  is  a good  approxi- 
mation for  the  range  of  interactions  < 3,  while  the  effective  mass 
m#  = mc(l  + p^w/2)  is  less  than  the  true  value. 

We  shall  study  the  problem  where  g(k)  = = constant  in  which  L is 

the  length  of  the  quantization  box.  For  the  one-dimensional  problem 


2(..) 

k 


L(.  . )dk 


If  one  introduces  dimensionless  units  by  measuring  energy  in  terms  of  hu> 
and  momentum  in  units  of  (2mtiu)^,  the  measure  of  coupling  strength  be- 
comes p^w  where 


.2  1 l 2/  h \/2mw\l 

p w ■ i TTTi  ry 

(hv) 


In  addition,  the  quantity 


2Mu 


lr. 


2hw  — -£ — (X  ♦ k + x)^ 
2m  ' ' 


has  the  dimensions  of  energy,  and  in  dimensionless  units  B = 


f2-E 


The  Fock  equations  including  three-quantum  terms  become 


( y 

VE  - x-y  X = P 
l Jo 
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|E-l-(\  + k)1 2j  ?i(k1)  = PX0+  ^2<klk2>  dkZ 

C2(klk2)=-A  [^(k,)*  5^)] 

r+oo  • 

♦ /3P  / £ (kjk^jJdk, 

J-oo  3 

|E—  3 — (X  ♦ kj  ♦ k2  ♦ k3)2J  C3(kjk2k3)  = ^^(kjkJ 


+ C (k2k3^  + £ ^k3k  1^1  * 

1 2 J 

► ' L ^2  f ► ' L h 

We  have  introduced  the  amplitudes  ^ = r-t*  s . s = -5—  b . etc. 

1(2*)^  1 2 Zn  2 

Let  us  now  discuss  these  equations  according  to  the  number  of  quanta.  The 


one -quantum  solution  is 


r ' P X 

C (M is — 

1 E - 1 - (X.  + k)‘ 


— E = -fifi- 
l/l  - E 


In  Fig.  2 we  plot  E(\  = 0)  = E as  a function  of  p it.  To  find  the  mean  velocity 

dE  t 

v we  must  compute  


1 + p w/2(l-  E)' 


t It  can  be  shown  by  methods  similar  to  those  in  Ref.  30  that  v = dE/dX. 
and  m#  = l/d2E/dX.2.  Compare,  however,  the  elegant  proof  of 
Bogoliubov  * 1 . 


Interaction  strength  P*  * 


Fig .2.  Energy  depression  vs.  interaction  strength. 
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To  find  the  effective  mass  for  small  values  of  the  total  momentum  we  set 
E = Eq  + (\^/m*)  and  expand  to  order  X2.  Then 


« = i + —Kl/L 
0 -EJ 


(a?) 


This  is  plotted  as  a function  of  P it  in  Fig.  3;  m#-^  1 + as  p w 0 and 
2 


m*  1.  50  as  p it 

To  find  the  two-quantum  solution  one  must  solve  the  pair  of  integral 
equations 

{E.,n  x *n  f*mt  wai. 

Jo  J-  oo  ( ( 18) 

Je  - i - (X  + k)2  + b)  C (k)  = p X ♦ P2  f (-I  — — . 

I J o J- oo  E-2-(X  + k + kr 

The  Liouville-Neumann  series  is.  with  B = and  the  abbreviations 

a = /l  - E - B , C = jZ  - E . 

r -ox.  f , rz 

a2  + (X  + k)2  c'-»o{a2  + + k')2J  {c2  + (k  + k + k')2J 


♦ P 


Jk 


j4  

wf-«o  |a2  + (X  + k')2|  |c2  + (X  + k ♦ k')2| 

r~ ♦ ...  i 

J- OO  ja2  + (X  + k')2j  |c2  + (X  + k + k')2j  j 

With  the  help  of  the  integral 


(•*9) 


a * ±) 


r *°°  dk' _ Ha  + T 

7-«o  (a2  + k'2)  (C2  + (k  + k’)2)  k2  + (a  +C) 


2 ’ 


(50) 


we  obtain 
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{ • — ~~~7~ — T f‘  * * ’Mi  * cK 

(X  + k)2  + a2  l k + (a  ♦ C)2  Va  C/ 

+ k + k’)2  ♦ C*  {(X  ♦ k*)2  + a2  . k'2  + (a  ♦ C)Z  ^ 


A i+c)(s4x*i)  . /a2_,Vl  . 1\2  1 

+ (p  w)  "xrMc72^2"  Mp  ^-q- 


(51) 


x2-e=0 


i + 2(2a  * c)(i4t:  f r)1 

(a  + c)2  [x2  + (2a  + C)2]  [x2  + (2a  + C)2]  2 


♦ . 


(52) 


In  order  to  evaluate  the  next  term  one  must  compute  the  three-denominator 
integral.  This  is  easy  to  do  but  will  not  be  written  here.  For  the  one- 
dimensional  problem  the  kernel  is 

1 


Kt*^  ) = 


1 


• <»•  ♦ * ♦£  I2  * C2}  [(X.f)2,a2] 


(53) 


4 2 3 3 

and  the  convergence  condition  is  § k € 4 C a . Proceeding  as  in  the  one- 
quantum  case,  we  find  for  the  effective  mass 


m*  = 1 + 


ZaC  ( 
o o 


* t + £)*  aifVcy  {-^r  ♦ (2  *fyi\ 

j o'  Loo  oo  o J 

^2  {*T  +i(2  + ^r)  * 2a  ♦ C V 

2a  + C ) I o C o o 1 

o o'  J 


(54) 


The  results  of  the  calculations  for  the  two-quantum  case  are  plotted  in  Figs. 

2 and  3.  It  is  seen  that  the  two-quantum  correction  to  the  energy  becomes 
2 

appreciable  for  {J  w > 1,  while  the  correction  to  the  effective  mass  is  im- 
portant for  (32w  > 0.3.  For  the  higher  values  of  p2w  it  is  necessary  to  in- 


• I • • 


• . • J 


^ • V . 


elude  several  higher  terms  in  the  series  (Eq.52).  For  example,  for  p it  = 2 
the  energy  is  Eq  = - 1.82  and  is  divided  as  — (1.51  + 0.25  + 0.05  0.01).  Even 

for  small  values  of  P^it  it  is  important  to  include  higher  terms  for  large  values 
of  X.  It  is  interesting  to  note  from  Eq.(48b)  that  the  onset  of  the  continuum  is 
depressed  to  an  energy  determined  by  1 — E — B = 0.  This  is  illustrated  in  Fig. 
4 where  we  plot  the  energy  vs.  momentum  for  = 1.  For  higher  X the  effec- 
tive mass  appears  to  be  large.  It  is  not  clear  how  higher  quanta  will  modify 
this  effect. 

In  order  to  estimate  the  energy  when  the  wave  function  includes  three 
quantum  contributions  it  is  necessary  to  eliminate  £ and  to  obtain  an  in- 
tegral equation  for  £ . One  finds 

2 


(e-2- 


(x  + kj  + k2r  ♦ a 


} Vkiki,  = A (Vk,)t  V“2,}+pi 


f £^k2k 3^  * £,(k3kl^ 

J E-  3 -l\  + k.  + k,  + 


where  A = 


f3  - E 


In  addition,  we  have  the  first  two  equations  of  set  (44).  We  can  no  longer 
find  a single  integral  equation  which  can  be  solved  by  iteration  and  whose 
convergence  can  be  demonstrated.  The  procedure  will  be  less  systematic. 
We  first  drop  the  integral  terms  above  and  find  for  £ (k.k?) 

1 ( 2 1 £ 

' { £<ki> + 5 (k2) 

£ = £ f r (56 

z {Z  J E - 2 - (X  + kj  + kzr  + A 1 

This  expression  differs  from  the  previous  one  only  by  the  presence  of  the 
constant  A.  The  wave  functions  and  energy  levels  can  therefore  be  found 
from  Eqs.(51)  and  (52)  by  replacing  a and  C by 


a = t/l  — E — D . C = *2  - E — A 


0.73 


Energy  vs.  momentum  for 
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where 


D = 


-2 

P * 

/2  — E — A 


I»  (k,k,)  found,  into 

The  next  correction  may  be  found  by  inserting  the  ^ 1 £ 

2 

the  integral  in  Eq.(55).  We  find 

t I 

• 2 /"+••  £ (k2k3^  * t (k3kl* 

6t(k,kJ= ^ i - 2 dk  . (58) 

' 1 1 - + k,r  ♦ C 4 j- oo  3 - E ♦ (k  + k.  + k3  + X)4 


x(kt  + k/ 


Since  C (k,M  can  be  expressed  explicitly  in  terms  of  kjk2  and  ^ this 

integral  can  be  evaluated.  Then  inserting  6 £ into  Eqs.  (44a)  and  (44b), 

2 

we  find  a correction  to  the  energy  of 


/. 


+ «© 


dk, 


oo  E - 1 - (X  + kj) 


/. 


dk. 


(X  + kj  + k£)2  + C'* 


^ (k2kj)  + ^ (k^kj) 

— 2 dk3 

3 - E ♦ (X  ♦ k ♦ k2  ♦ k3)‘ 


(59) 


where 


I I 

( C O',)  * £ O'  ) > 

1 1 1 

2 - E + (X  + kj  + k2)2  - A 


and  £ may  be  obtained  from  Eq.(51).  Higher  order  corrections  are  ob- 
1 

tained  by  iteration.  From  Figs.  2 and  3 we  see  that  three -quantum  effects 
influence  the  energy  and  effective  mass  to  an  important  extent  for  p^w  > 3. 
The  effective  mass  is  particularly  sensitive  to  the  number  of  quanta  included 
and  for  p^w  > 3 it  is  possible  that  the  method  of  including  a limited  number  of 
quanta  breaks  down.  The  calculations  of  Fig.  3 seem  to  point  to  a breakdown 
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of  the  linear  dependence  of  m#  on  p2w  for  p2ir  > 3. 


5.  Three-Dimensional  Case 


Let  us  now  examine  the  actual  three-dimensional  case.  We  introduce 
dimensionless  units  by  measuring  energy  in  units  of  hu>  and  momentum  in 
units  of  j^zinhu  so  that  the  relation  E =■  \2/2m  for  a free  particle  goes  over 


to  E = X2.  We  make  use  of  the  fact  that  £ ♦ — f k2dkdX2.  Then 

k (2n)3  J 

with 

...  . . e2  /2mu^2 

a 

l \ . 

n 


1/1  1 \ e2  ( Zmu>\] 

= 2 V~T"  7)lZ>\-T-) 


the  equation  for  the  energy,  including  two  quantum  terms,  becomes 


{X2_E}=I7  Iff  k)*t?E-B} 

dk'd  £1  * 


1 + 


Zl T 


/// 


dk'd  afc 

j (L  + k + Is')2  + 2 - e|  |(X  + k')2  + 1 - E - B j + • 

* 


where 


B (&•  *)  = fx,  V~1T]  ain  l (j — * kl; — )• 

*-  -■  Kiz  + (X  + k)2  - E ' 

Neglecting  the  two-quantum  contributions  one  finds  the  one-quantum  ex- 
pression 

X 


(60) 


(61) 


.2  _ a . -1 

X — E = y sin 


/ 1 + X2  — E 


(62) 


The  lowest  energy  level  is  obtained  for  X.  = 0 when  — Eq  = a// 1 — Eq  . This 
is  the  same  expression  that  we  had  in  the  one-dimensional  case  provided  v>e 
identify  the  p2ir  used  there  with  the  a used  here.  One  can  ob-  mi  the  effective 
mass  by  writing  E = Eq  + X.2/ m#>  and  developing  in  powers  of  X.  W «_  fir.-* 


' ■*  . w 
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m 


1 + 


(i  - 


E ) 
o' 


(63) 


The  energy  Eq  is  plotted  as  a (unction  of  a in  Fig.  3.  In  the  limit  of  vanishing 
a one  finds  m#  1 ♦ a/6  so  that  the  correction  is  l/J  that  for  the  one- 
dimensional case.  For  very  large  a,  Eq  “►  o^  and  m# -►  1.13. 

Let  us  now  turn  our  attention  to  the  two-quantum  contributions.  The 
results  of  the  one -dimensional  case  indicate  that  for  actual  crystals  a two- 
quantum  solution  will  not  be  adequate.  Since  the  calculations  are  involved, 
we  shall  restrict  ourselves  to  a discussion  of  the  energy  for  X.  = 0 as  a 
function  of  interaction  strength.  The  results  for  the  energy  are  similar  to 
those  for  the  one-dimensional  case.  For  the  effective  mass  we  would  expect 
approximate  agreement  with  the  law  m*=  m ( 1 + a/6). 

We  notice  that,  because  of  the  presence  of  the  term  B(k)  in  the  de- 
nominators, one  does  not  have  straightforward  integrals.  Consider  the  de- 
nominator k^  * 1 — E — B with  B = a/k  sin  * k/J^2  + k^  — E.  For  k = 0, 

B — ► a// 2 — E , while  for  large  k,  B 0.  It  can  be  shown  that  assuming 
B = const.  = a// 2 — E introduces  errors  of  the  order  1 percent  in  the  energy. 
In  addition,  a better  constant  value  for  B may  be  determined  after  one  has  ap- 
proximate values  for  E.  Thus  the  first  term  in  the  Liouville -Neumann  series 
gives 

X^  — E = sin  * - - — (64) 

^1  + X2  - E — B 


In  order  to  compute  the  second  term  in  the  series,  one  must  evaluate 
an  integral  of  type 


dk'd 


r7^ 


I 


(65) 


•*  • ' 


...... 


<"■  < * • V . 


where  we  have  put  = 1 — E — B,  b ^ = 2 — E . We  write  I in  the  form 


-IlUIf  ^ 


dk  dk 


(k*  + a)  |(k  - k’)Z  + t>  l (k,Z  ♦ a1)  kV' 


and  consider  first 

J -Jtf : T 5 = /8(k';  Mi  - i')  dk\ 

JJJ  k'2(k'2  ta2)  ((k  -k')2  * b2)  ^ 


where 


e^')  • . ,2.. .z ‘ z.  '•  Mi-i')’.  rr. 


k'‘(k'‘  t a ) 


We  now  make  use  of  the  theorem 


(k-k,’)4-  ♦ b* 


/©  (P)  <t>(P)*~'P  ' * dP  = / 9(k')  *(k-k’)  dk'  . 


where 


is  the  Fourier  transform  of  9(k  ) and  $ip)*fn  c'bp/p  ■»  the  Fourier 
transform  of  $(k').  We  have 

2 -III 

Introducing  the  angle  cos  9 between  k and  p , integrating  over  9,  and 
inserting  in  the  expression  for  I,  one  finds 

T - r dk  /"*,  „ ~X>  b ( 1 e ' Pa) 

1 - -r  J zrr~r  / dp*H  - — - — is‘nk p , 

a k + a */0  p k 


and  finally, 


3 


An  expression  for  I can  be  obtained  in  the  form  of  a series  by  writing 


Fig.  5.  Energy  of  lowest  state  (three  dimensions) 


» J ' • 
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. -1 
tan  x 


x — 


13.15 


2 


x 


< 1 


X 


2 


> 1 . 


Thus,  to  the  second  term,  we  find  (for  k = 0)  — E = a //l  — E — B + I(a,  b). 
The  contributions  of  higher  terms  in  the  Liouville-Neumann  series  become 
progressively  more  difficult  to  evaluate.  However,  each  term  gives  a 
positive  contribution  to  the  right-hand  side  of  Eq.(60).  Furthermore,  one 
can  over-estimate  the  contribution  of  each  term  by  replacing  the  coupling 
denominator  (k  + k ) I b by  b . Then  each  term  is  less  than 


Q 1 

(2-  E)  /l  - E - B 

times  the  previous  one,  and  if  one  uses  the  overestimate,  one  obtains  values 
of  E which  are  too  low.  On  the  other  hand,  neglecting  higher  terms  entirely 
gives  values  of  E which  are  too  high. 

The  results  of  the  calculations  of  E vs.  a are  presented  in  Fig.  5. 

The  behavior  is  very  similar  to  that  in  one  dimension. 


